% 2021.10.11
% 本程序是HW1改进的MIMO程序

% 参数：1200kg、(200, 200, 150)kN/m、0.02, 0.02, 0.04
% /1000 *1000不同？
m = 1200;
mass_m = m * eye(3);
stif_m = 1000*[400 -200 0;
    -200 350 -150;
    0 -150 150];
% mass_m\stif_m
% 阻尼比
cri_damp_m = [0.02, 0.02, 0.04];
% 自振频率
% inv(mass_m)*mass_m=eye(3)
% [eig_vec, eig_val] = eig(inv(mass_m)*stif_m);
[eig_vec, eig_val] = eig(mass_m\stif_m);
% 如何对应频率
omega_m = sort(sqrt(diag(eig_val)))';     % 5.6419   14.7706   22.3607
